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' 1. Introduction 

oo : 

O I It is now accepted that QCD is the theory which underlies all of nuclear structure. A 
Q ' fundamental question we may then ask is: How do nuclear energy levels change as we 
■ vary the quark masses in the QCD Lagrangian? This is not solely an intellectual exercise; 
^ . there are hints that the fundamental parameters of the standard model, such as and 
^ ', nig, may be time-dependent [Hj. The successful predictions of Big-Bang nucleosynthesis 
r-^ I (BBN) can be used to constrain the time dependence of these parameters and thereby 
^ ' search for physics beyond the standard model [12] [El- However in order to do so one must 
C ■ know how nuclear physics depends on the fundamental parameters. A more practical 
>• . motivation for understanding the quark mass dependence of nuclear physics, and more 
^ ! generally of hadronic physics, is that in the near future it is through lattice QCD simula- 
^ [ tions that definitive predictions in nuclear physics will be made directly from QCD. And 
- - ' because lattice QCD is currently simulated with unphysically large quark masses, there is 
an inevitable extrapolation to physical quark masses. If rigorous predictions are to be had 
from the lattice, this extrapolation must be controlled in a precise way. Fortunately, the 
quark-mass dependence of few-nucleon systems can be studied by exploiting hierarchies 
of scales: i.e. by using effective field theory (EFT) [ E] [ [ El • Given the current state 
of technology, simulations of multi-nucleon systems are intractable, but realistic simula- 
tions of two-nucleon systems are feasible. Arguably, the most promising method is to 
calculate scattering phase shifts directly using Liischer's finite-volume algorithm, which, 
for instance, expresses the ground-state energy of a two-particle state as a perturbative 
expansion in the scattering length divided by the size of the box [ 7]. This method has 
been used successfully to study hit scattering [|B]. Attempts have been made to compute 
nucleon-nucleon (NN) scattering parameters in lattice QCD; Ref. [ Q\ computes the ^5*0 
and ^5*1 scattering lengths in quenched QCD (QQCD) using Liischer's method. Of course 
in the NN system, the scattering lengths are much larger than the characteristic physical 
length scale given by Ag^^,. Nevertheless one may expect that at some unphysical value 
of the pion mass used in a lattice simulation, the scattering lengths relax to natural val- 
ues, thus allowing their determination from the lattice. This is yet another motivation 
for understanding the quark-mass dependence of few-nucleon systems. 
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Presently, unquenched lattice simulations with the physical values of the light-quark 
masses are prohibitively time-consuming, even on the fastest machines. While some 
quenched calculations can be performed with the physical quark masses there is no limit in 
which they reproduce QCD, and consequently they should be considered to be a warm-up 
exercise for the "real thing". Relatively recently it was realized that partially-quenched 
(PQ) simulations, in which the sea quarks are more massive than the valence quarks, 
provide a rigorous way to determine QCD observables and are less time-consuming than 
their QCD counterparts. It is with PQ simulations that nuclear-physics observables will 
first be calculated rigorously from QCD. As PQQCD reproduces QCD only in the limit 
in which the sea-quark masses become equal to the valence-quark masses which, in turn, 
are set equal to the physical values of the light quark masses (we call this the QCD limit), 
there are some interesting features of the PQ theory that are distinct from nature away 
from the QCD limit. In QCD, the long-distance component of the NN potential is due 
to OPE, as discussed above. However, in PQQCD there is also a contribution from the 
exchange of the ?7-meson (in the theory with two flavors of light quarks). In QCD such 
an exchange is suppressed due to the large mass of the rj compared to the vr. However, in 
PQQCD the rj propagator has a double-pole component that depends on the pion mass 
due to the hairpin interactions with a coefficient that depends upon the difference be- 
tween the masses of the sea and valence quarks. Therefore, away from the QCD limit the 
long-distance component of the NN potential is dominated by one-eta exchange (OEE) 
and falls exponentially with a range determined by the pion mass, ~ e"™'"'', as opposed 
to the familiar Yukawa type behavior [KlUj. All is not lost however: it is straightforward 
to develop the partially-quenched EFT which matches to a partially-quenched lattice 
simulation [ITT]. 

A second approach to the NN system on the lattice is to study the simplified problem 
of two interacting heavy-light particles [112]. It has been suggested that lattice QCD 
simulations of the potential between hadrons containing a heavy quark will provide insight 
into the nature of the intermediate-range force between two nucleons [IE]- While the NN 
potential is not itself an observable, one may instead consider heavy systems. In the 
heavy-quark limit, the kinetic energy of the heavy hadrons is absent and the lowest-lying 
energy eigenvalues, which can be measured on the lattice, are given by the interaction 
potential. All discussions to date have addressed two heavy mesons. This case is somewhat 
complicated by the fact that there are degeneracies in the heavy-quark limit which require 
a coupled- channel analysis. By contrast, the AqAq interaction (where the Ag is the lowest- 
lying baryon containing one heavy quark, Q), does not suffer from this complication [ I14j . 
Moreover, since the Aq is an isosinglet, there is no OPE, and the leading large-distance 
behavior is governed by two-pion exchange (TPE), which is physics analogous to the 
intermediate-range attraction in the NN potential. 

In this proceedings I will review recent work in developing the EFTs relevant to under- 
standing the quark-mass dependence of the NN system and the AqAq potential, both in 
QCD and in PQQCD. 
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2. The Quark Mass Dependence of NN 

During the last decade significant effort has been expended in constructing an EFT 
to describe nuclear physics. While it is straightforward to write down all possible terms 
in the effective Lagrangian for two or more nucleons, arriving at a consistent power- 
counting has proved to be a difficult task. Weinberg's (W) original proposal [ il5j for 
an EFT describing multi-nucleon systems was to determine the NN potentials using the 
organizational principles of the well-established EFTs describing the meson-sector and 
single-nucleon sector (chiral perturbation theory, xPT), and then to insert these potentials 
into the Schrodinger equation to solve for NN wavef unctions. Observables are computed 
as matrix elements of operators between these wavef unctions. W power-counting has 
been extensively and successfully developed during the past decade to study processes 
in the few-nucleon systems [ J^. This method is intrinsically numerical and is similar 
in spirit to traditional nuclear-physics potential theory. Unfortunately, there are formal 
inconsistencies in W power-counting [^\, in particular, divergences that arise at leading 
order (LO) in the chiral expansion cannot be absorbed by the LO operators. Problems 
persist at all orders in the chiral expansion, and the correspondence between divergences 
and counterterms appears to be lost, leading to uncontrolled errors in the predictions for 
observables. This formal issue was partially resolved by Kaplan, Savage and Wise (KSW) 
who introduced a power-counting in which pions are treated perturbatively [ ^] . The 
NN phase-shifts and mixing angle in the ^5*0 and ^Si — ^Di coupled-channels have been 
computed to next-to-next-to-leading order (N^LO) in the KSW expansion by Fleming, 
Mehen and Stewart (FMS) [E] from which it can be concluded that the KSW expansion 
converges slowly in the ^5*0 channel and does not converge in the ^5*1 — ^Di coupled- 
channels. Therefore, neither W or KSW power-counting provide a complete description 
of nuclear interactions. 
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Figure 1. The left (right) panel shows 
channel) function of the pion mass, 
and the black region corresponds to r] = 
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the scattering length in the ^S'o-channel (^5*1- 
The light gray region corresponds to 77 = 1/5 
1/15. In the ^S'l-channel the parameter diQ is 
: dw < -0.17 GeV-2 and rfis = -0.51 GeV^l 



The problems with W and KSW power-counting appear to have been resolved in Ref. [ 
1^, which from this point on I will refer to as BBSvK. It was realized in FMS that the con- 
tributions to the amplitude that lead to non-convergence in the ^Si — ^Di coupled-channels 
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persist in the chiral limit (it is the chiral hmit of iterated one-pion-exchange (OPE) that is 
troublesome). Therefore, in BBSvK power-counting the scattering amplitude is an expan- 
sion about the chiral limit. This recovers KSW power-counting in the ^5*0 channel, where 
FMS found it to be slowly converging. However, in the ^5*1 — '^Di coupled-channels, the 
chiral limit has contributions from both local four-nucleon operators and from the chiral 
limit of OPE. It is these two contributions that must be resummed using the Schrodinger 
equation to provide the LO scattering amplitude in the ^S'l — ^Di coupled-channels. 

In recent papers by Savage and the author [ ^ and also by Epelbaum, Glockle and 
MeiBner [ [Bj EFT was used to determine the mg-dependence of scattering in the two- 
nucleon sector. Remarkably, in the ^S'o-channel KSW power-counting can be used to 
derive an analytic expression for the scattering length, 

where 7 is a LO constant and D2{fi) is a combination of coefficients of operators with a 
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Figure 2. The deuteron binding energy as a function of the pion mass (from the phys- 
ical value to the chiral limit). We do not show values below m^r = 60 MeV since the 
deuteron can be both bound and unbound. The shaded region corresponds to 77 = 1/5, 
-2.61 GeV-^ < rfi6 < -0.17 GeV"^ and -1.54 GeV'^ < Sig < -0.51 GeV'^ 

single insertion of rriq. Unfortunately, D2 is presently unknown (in the ^5*1 — ^Di channel 
D2 contributes to, for instance, vr-deuteron scattering, however only at one order beyond 
the current state-of-the-art [EDlfE])- The best that one can do at present is to use naive 
dimensional analysis (NDA) to estimate a range of reasonable values for D2, defined by 
a parameter -C 1 [E]. The results of NDA are shown in Fig. [H We use scatter plots as 
the point density represents the probability associated with a particular set of low-energy 
constants. NDA suggests that the di- neutron remains unbound in the chiral limit, while 
a relatively small increase in rriq could lead to a bound di-neutron. 

In the ^5*1 — ^Di coupled channels the situation is somewhat more complicated. At 
NLO in BBSvK counting there is not only OPE, but also the chiral limit of TPE. As a 
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consequence, there are additional counterterms in the single nucleon sector that contribute 
in this channel but do not contribute in the ^5*0 channel, in particular dig and diQ associated 
with the pion-nucleon interaction, and I4 associated with This is in addition to the 
D2 contribution in the ^5*1 channel. The allowed regions for dis and rfig are given in Ref. [ 
122] J and I4 is known. Fig. ^ shows the presently allowed values of the scattering length 
in the ^5*1 channel where we again have used NDA to estimate the possible values for D2. 
It is clear that for the range of parameters considered the deuteron could be bound or 
unbound in the chiral limit, and at present one cannot make a more definitive statement. 
This last statement is in disagreement with the conclusion of Ref. [EI; a discussion of this 
disagreement is given in Ref. The quark- mass dependence of the deuteron binding 
energy is shown in Fig. |21 Notice that the ^Si scattering length relaxes to natural values 
of ~ 1 fm as the pion mass is increased beyond ~ 200 MeV (see Fig. P). One anticipates 
similar behavior in the partially-quenched theory. Given current uncertainties in strong 
interaction parameters, particularly D2, it is at present unclear whether the same is true 
in the ^5*0 channel (see Fig. H]). 

As an interesting application of these results, Ref. [ |2] derive constraints on the time 
variation of the Higgs vacuum expectation value {(p) through the effects on BBN, including 
the effect of the change in the deuteron binding energy, which alters both the ^He and 
deuterium abundances significantly. See Fig. El 
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Figure 3. The solid curves represent the constraints from the ^He and deuterium abun- 
dances, assuming no change in the deuteron binding energy (left panel) and a change in 
the deuteron binding energy as shown in Fig. |21 The region allowed by BBN is shaded. 



It is quite simple to construct the partially-quenched effective field theory from the 
known QCD results. And too, it is gratifying to see that one can obtain analytic results 
for many observables in the ^5*0 channel and in the higher partial waves. For instance, the 
^Pi scattering volume is given by 

n('P)= ^^^^ I ^"^^ ^Is-^l 



where go is an axial coupling and mss is a Goldstone boson mass containing two sea 
quarks. Notice that the QCD limit agrees with the well-known results [123- One should 
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be concerned about the range of sea and valence quark masses for which this theory 
converges. In QCD it is found that the NN EFT converges for and momenta less than 
of order A^vat ~ 300 MeV, and one suspects that the same radius of convergence will exist 
in the partially-quenched theory. If this is indeed the case, lattice calculations will be 
required with meson masses of less than ~ 300 MeV in order to match to the EFT and 
use it to make predictions about nature. This is somewhat more restrictive than in the 
meson and single nucleon sectors and therefore one would like to see convergent results 
in those sectors before being confident in results obtained in the multi-nucleon sectors. 

3. The AqAq Potential 

The lowest-lying baryons containing a single heavy quark can be classified by the spin of 
their light degrees of freedom (dof), si, in the heavy-quark limit, mq — >■ oo. Working with 
two light flavors, u and d quarks, the light dof of the isosinglet baryon, Aq, have si = 0, 
while the light dof of the isotriplet baryons, T,q^'^ and Sq^'°*, (the superscript denotes 
the third component of isospin) have s/ = 1 ^. In the heavy-quark limit the spin-i 
baryons are degenerate with the spin-| Sq^'°* baryons, but are split in mass from the Aq 
by an amount that does not vanish in the chiral limit. As the light dof in the Aq have s/ = 
in the heavy-quark limit, the light-quark axial current matrix element vanishes, and thus 
there is no AqAqtt interaction at leading order in the heavy quark expansion. This means 
that there is no OPE (or OEE) contribution to the AqAq potential in QCD and PQQCD, 
and therefore there is no long-distance component in the AqAq potential. It is the TPE 
box and crossed-box diagrams that provide the longest- distance interaction between two 
Aq's. In addition, there are local four-Ag operators at the same order in the chiral 
expansion but such local interactions give coordinate-space delta-functions. Analogous 
diagrams in the two-nucleon sector provide part of the intermediate-range component of 
the NN potential. However, it is important to realize that there are additional interactions 
that contribute to the intermediate-range component of the NN potential in the chiral 
expansion, for instance contributions from the Weinberg- Tomazawa term and also from 
higher- dimension NNtttt vertices. Therefore, while the AqAq potential provides a window 
into the nature of the intermediate-range NN interaction, it certainly does not provide a 
complete description. 

If the Aq and Sg^ were degenerate, we would be required to solve the coupled-channel 

system with AgAg and Sg-'Sg^ coupled to / = 0. In the charmed sector the Sc — Ac 
mass splitting is A = 167.1 MeV and the S* — Ac mass splitting is A = 232.7 MeV, 
and we use the spin- weighted average of these splittings to estimate A ~ 211 MeV in 
the heavy-quark limit. There is no symmetry reason for this mass-splitting to vanish in 
the chiral limit, and hence there is no infrared divergence that requires a coupled-channel 
analysis. In the power-counting we treat A ~ 771,^ and take the Maq , A^Sq — ^ oo limit 
in evaluating the diagrams in Fig. |3 With this power-counting, one can directly use the 
Feynman rules of Heavy-Baryon Chiral Perturbation Theory (HBxPT) to describe the 
low-momentum dynamics of the nucleon and A-resonance, without the need to resum the 
baryon kinetic energy term as is the case for the box and crossed-box diagrams in the 

^For three light flavors the baryons fall into a 6 ® 3 of SU(3). 
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nucleoli sector. Evaluating the diagrams in Fig. |3 and then Fourier transforming them to 
position space is straightforward [E]- For asymptotically large distances, the potential 
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Figure 4. The box and crossed-box diagrams that give the longest-distance component 
of the AqAq potential. Heavy-quark symmetry forbids AqAq intermediate states in the 
box and crossed-box diagrams. 

is well-approximated by 

o 4 9/2 

16 7r5/2 r5/2 A /4 

which exhibits the expected fall off with a length scale set by twice the mass of the pion 
and where the dots represent subleading contributions in the large-distance expansion. 

The extension of the heavy-baryon sector from QCD to PQQCD is straightforward 
and the leading modifications to the AqAq potential are easily computed [ [T3] . Fig. El 
exhibits the potential as a function of m^y, the mass of the Goldstone boson consisting of 
a valence quark and a sea quark, for two values of r. One may wonder whether "hairpin" 




mj.y [MeVJ mj.y [MeV] 

Figure 5. The left panel shows V^'^{r) evaluated at r = 1 fm as a function of the meson 
mass msV) while the right panels shows V^^{r) evaluated at r = 2 fm. The vertical axis 
is in units of MeV. When msv = ''^Tr the value of V^'^{r) is equal to V'^'^^ir). 



contributions may enter in such a way as to dominate the potential. The leading hairpin 
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Figure 6. The leading hairpin contribution to the AqAq potential. 



contribution enters through the diagram in Fig. IHl At asymptotically large distances this 
contribution to the potential becomes 

^ ^"^^ 64 7r5/2 A6 ^ + • • • ■ 

where the dots represent subleading contributions in the large-distance expansion. While 
at asymptotically large distances this contribution is larger than that from the box and 
crossed-box diagrams, asymptopia finally sets in at distances at which all contributions 
are numerically insignificant. 



4. Conclusions 

Understanding how nuclei and nuclear interactions depend upon the light-quark masses 
is a fundamental aspect of strong-interaction physics. Recent work has been able to 
explore the m^-dependence of two-nucleon systems using a recently-developed effective 
field theory and naive dimensional analysis. In the ^^o-channel we expect that di-nucleon 
systems, such as the di-neutron, are unbound for all values of ntq less than their physical 
values. However, for larger than their physical values both bound and unbound 
systems are presently consistent with data and NDA. In the ^Si — ^Di coupled-channels, 
where the deuteron resides for the physical values of the quark masses, the deuteron 
may or may not be bound in the chiral limit. A more definitive statement can only be 
made with a more precise determination of the nN coupling diQ and a determination of the 
coefficients of the leading rrig-dependent four-nucleon operators, D2. As discussed in Ref. [ 
it is likely that a determination of D2 will require a future lattice QCD calculation. 

Recent work has computed the potential between two Ag's at leading order in effective 
field theory in both QCD and PQQCD. The size of the leading contribution from hairpin 
interactions in PQQCD has been estimated. Evidently the partially-quenched AqAq 
potential does not suffer from some of the (partial-) quenching problems that plague the 
NN potential due to the absence of single pseudo-Goldstone exchange. The computed 
potentials will allow for the chiral extrapolation of lattice calculations performed with 
unphysically large sea quark masses. As these potentials fall off with a mass scale set by 
~ 2m,r, they are quite small for baryon separations greater than r ~ 1.5 fm. Therefore, 
the theoretical advantages of studying this system to learn about the NN potential may be 
undermined by the difficulties in extracting a signal from lattice simulations. However, 
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the simplifications introduced by only having two light quarks, and a single infinitely- 
massive quark to fix the inter-baryon separation makes this system a prime candidate for 
studying inter-baryon interactions. It is very exciting indeed to be so close to making 
fundamental statements about nuclear physics. 

I would like to thank the organizers for a very enjoyable meeting and my collaborators 
Daniel Arndt and Martin Savage for sharing their insight. This work is supported in part 
by the U.S. Department of Energy under Grant No. DE-FG03-00-ER-41132. 
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